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Abstract — We  investigate  spatial  interference  statistics  for 
multiGigabit  outdoor  mesh  networks  operating  in  the  unlicensed 
60  GHz  “millimeter  (mm)  wave”  band.  The  links  in  such 
networks  are  highly  directional:  because  of  the  small  carrier 
wavelength  (an  order  of  magnitude  smaller  than  those  for  existing 
cellular  and  wireless  local  area  networks),  narrow  beams  are  es¬ 
sential  for  overcoming  higher  path  loss,  and  can  be  implemented 
using  compact  electronically  steerable  antenna  arrays.  Direction¬ 
ality  drastically  reduces  interference,  but  it  also  leads  to  “deaf¬ 
ness,”  making  implicit  coordination  using  carrier  sense  infeasible. 
In  this  paper,  we  make  a  quantitative  case  for  rethinking  medium 
access  control  (MAC)  design  in  such  settings:  unlike  existing 
MAC  protocols  for  omnidirectional  networks,  where  the  focus  is 
on  interference  management,  we  contend  that  MAC  design  for 
60  GHz  mesh  networks  can  essentially  ignore  interference,  and 
must  focus  instead  on  the  challenge  of  scheduling  half-duplex 
transmissions  with  deaf  neighbors.  Our  main  contribution  is 
an  analytical  framework  for  estimating  the  collision  probability 
in  such  networks,  as  a  function  of  the  antenna  patterns  and 
the  density  of  simultaneously  transmitting  nodes.  The  numerical 
results  from  our  interference  analysis  show  that  highly  directional 
links  can  indeed  be  modeled  as  pseudowired,  in  that  the  collision 
probability  is  small  even  with  a  significant  density  of  transmitters. 
Furthermore,  simulation  of  a  rudimentary  directional  slotted 
Aloha  protocol  shows  that  packet  losses  due  to  failed  coordination 
are  an  order  of  magnitude  higher  than  those  due  to  collisions, 
confirming  our  analytical  results  and  highlighting  the  need  for 
more  sophisticated  coordination  mechanisms. 

Index  Terms — 60  GHz  networks,  millimeter  (mm)  wave  net¬ 
works,  interference  analysis,  medium  access  control,  wireless 
mesh  networks. 


I.  Introduction 

The  several  gigahertz  of  unlicensed  spectrum  available 
worldwide  in  the  60  GHz  “millimeter  (mm)  wave”  band, 
coupled  with  the  recent  progress  in  low-cost  mm-wave  radio¬ 
frequency  integrated  circuit  (RFIC)  design,  indicates  that  the 
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mass  market  deployment  of  multiGigabit  wireless  networks  is 
on  the  cusp  of  feasibility.  There  is  significant  industry  interest 
in  indoor  applications  of  this  technology  to  wireless  personal 
area  networks  [1] — [3],  wireless  local  area  networks  [4],  and 
wireless  uncompressed  HDTV  [5].  However,  use  of  the  60 
GHz  band  is  also  very  attractive  for  outdoor  mesh  networks 
with  relatively  short  link  ranges  of  100-200  meters:  oxygen 
absorption  in  the  60  GHz  band  produces  propagation  losses  of 
10-16  dB/km,  so  that  it  only  adds  about  3  dB  to  the  link  budget 
for  a  200  m  link.  Such  multiGigabit  outdoor  mesh  networks 
provide  an  easily  deployable  broadband  infrastructure  that  can 
have  a  multitude  of  applications,  including  wireless  backhaul 
for  picocellular  networks,  and  as  an  alternative  to  fiber  to  the 
home.  Our  goal  in  this  paper,  however,  is  to  point  out  that,  in 
order  to  realize  the  potential  for  such  networks,  it  is  necessary 
to  rethink  the  design  of  network  protocols,  accounting  for 
the  unique  physical  characteristics  of  mm-wave  links.  In 
particular,  we  observe  that  links  in  such  networks  must  be 
highly  directional,  and  provide  a  framework  for  analyzing  the 
spatial  interference  statistics  in  order  to  obtain  quantitative 
guidelines  for  medium  access  control  (MAC)  design. 

The  fundamental  distinguishing  feature  of  the  60  GHz 
band  relative  to,  say,  the  2.4  GHz  WiFi  band,  is  the  order 
of  magnitude  difference  in  wavelength.  For  omnidirectional 
communication,  free  space  propagation  loss  scales  as  A2, 
where  A  =  c/fc  is  the  carrier  wavelength,  with  c  denoting  the 
speed  of  light  and  fc  the  carrier  wavelength.  The  wavelength 
at  60  GHz  is  5  mm,  while  the  wavelength  at  2.4  GHz 
(WiFi  band)  is  12.5  cm,  so  that  the  propagation  loss  for 
omnidirectional  transmission  and  reception  is  625  times,  or 
28  decibels  (dB),  worse  at  60  GHz  than  at  2.4  GHz.  On  the 
other  hand,  for  a  given  antenna  aperture,  directivity  scales  as 
A.  Thus,  fixing  the  antenna  aperture  at  each  end,  we  gain 
by  a  factor  of  A,  making  the  overall  propagation  coefficient 
scale  as  i ;  this  corresponds  to  a  net  gain  of  28  dB  in 
going  from  2.4  GHz  to  60  GHz.  Thus,  given  the  difficulty  of 
producing  large  amounts  of  radio  frequency  (RF)  power  at 
mm  wave  frequencies  using  low-cost  silicon  implementations, 
employing  highly  directive  antennas  at  both  transmitter  and 
receiver  is  essential  for  the  range/rate  combinations  we  wish 
to  achieve.  Fortunately,  the  small  carrier  wavelength  also 
allows  for  the  realization  of  compact,  electronically  steerable 
antenna  arrays  (with  elements  realized  as  patterns  of  metal  on 
circuit  board),  as  demonstrated  in  [6],  [7],  While  hardware  is 
not  the  focus  of  this  paper,  we  note  that  it  is  possible  to  design 
nodes  which  provide  steerable  narrow  beams,  while  providing 
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TOP  VIEW 


(a)  Yagi-Uda  antenna  array 


(b)  Open-slot  antenna  array 


Fig.  1.  Steerable  arrays  with  high  directivity  are  realizable  with  compact 
form  factors.  Each  element  in  the  array  can  have  significant  directivity  (e.g., 
up  to  20  dBi),  and  combining  several  elements  with  overlapping  fields  of 
view  can  achieve  directivities  on  the  order  of  25—30  dBi. 


a  360°  field  of  view;  see  Fig.  1  for  a  concept  design.  Such 
omni-coverage  yet  highly  directional  nodes  are  the  analogue 
of  omnidirectional  nodes  at  lower  carrier  frequencies,  and 
enable  plug-and-play  deployment,  since  they  can  steer  highly 
directional  beams  to  neighbors  in  any  direction. 

Wireless  networks  at  lower  carrier  frequencies  typically 
employ  omnidirectional  (or  mildly  directional)  links,  so  that 
interference  is  the  main  performance  bottleneck.  Thus,  conven¬ 
tional  MAC  design  for  such  networks  focuses  on  mechanisms 
for  interference  management.  In  this  paper,  we  make  a  quan¬ 
titative  argument  that  a  drastic  rethinking  of  this  approach  is 
needed  for  the  design  of  MAC  protocols  for  60  GHz  mesh 
networks.  The  high  directivity  of  60  GHz  links  implies  that 
the  interference  among  simultaneously  active  links  can  be 
expected  to  be  significantly  smaller  than  for  omnidirectional 
links.  It  also  leads  to  deafness:  we  can  no  longer  count 
on  nodes  being  able  to  monitor  transmissions  destined  for 
other  nodes,  so  that  the  use  of  coordination  mechanisms 
such  as  carrier  sense  or  broadcast  control  messages  becomes 
infeasible.  We  argue,  therefore,  that  MAC  design  for  highly 
directional  60  GHz  networks  can  ignore  interference  up  to 
first  order  (e.g.,  by  being  reactive  rather  than  proactive  about 
interference  management),  and  should  focus  attention  instead 
on  the  problem  of  transmitter-receiver  coordination  in  the  face 
of  deafness.  While  our  focus  in  this  paper  is  to  quantitatively 
justify  these  new  design  guidelines  rather  than  to  present  a 
candidate  MAC  protocol,  we  note  that  preliminary  results  for 
a  distributed  MAC  protocol  designed  based  on  these  guidelines 
are  promising  [8], 

The  contributions  of  this  paper  can  be  summarized  as 
follows: 

•  We  analyze  the  statistics  of  the  spatial  interference  due 
to  multiple  simultaneous  and  uncoordinated  transmissions, 
with  transmitters  distributed  uniformly  on  the  plane  and  link 
orientations  chosen  randomly.  The  collision  probability  at  a 
designated  receiver  is  computed  as  a  function  of  the  node 
density  and  antenna  pattern,  and  the  analysis  accounts  for 
square  law  line-of-sight  propagation  loss  as  well  as  the  ex¬ 
ponential  loss  due  to  oxygen  absorption.  We  consider  both  a 
“protocol  model”  (where  a  collision  occurs  if  some  interferer 
exceeds  a  threshold)  and  a  “physical  model”  (in  which  the 


sum  interference  exceeds  a  threshold),  where  the  terminology 
is  borrowed  from  the  work  of  Gupta  and  Kumar  [9], 

•  Our  numerical  results  show  that,  for  the  antenna  directivi¬ 
ties  that  we  envision,  the  collision  probability  is  small  even 
without  any  coordination  among  transmitters.  This  motivates  a 
pseudowired  abstraction  to  serve  as  a  first-order  approximation 
of  mm-wave  wireless  links  for  the  purpose  of  MAC  design. 
Under  this  abstraction,  transmissions  on  different  links  do 
not  interfere  with  each  other.  However,  unlike  a  truly  wired 
node,  mm-wave  network  nodes  have  a  half-duplex  constraint, 
i.e.,  they  can  either  send  or  receive  at  a  given  time,  but  not 
both.  This  allows  MAC  designers  to  concentrate  on  developing 
lightweight  protocols  to  schedule  transmissions  in  a  “deaf” 
network  subject  to  half-duplex  constraints. 

•  We  verify  that  transmitter-receiver  coordination  is  indeed  the 
bottleneck  by  simulating  a  directional  slotted  Aloha  protocol 
for  60  GHz  mesh  networks  with  random  topologies.  The  simu¬ 
lation  results  show  that  packet  losses  due  to  failed  coordination 
are  an  order  of  magnitude  higher  than  those  due  to  collision. 

The  results  in  this  paper  build  on  the  preliminary  results 
presented  in  our  conference  publication  [10].  While  our  con¬ 
ference  publication  employed  Monte-Carlo  simulations  for 
estimating  collision  probabilities  for  the  physical  model,  a 
major  new  contribution  in  this  paper  is  that  we  provide  new 
analytical  approaches  that  yield  good  approximations  as  well 
as  tight  upper  bounds.  We  also  provide  a  more  extensive  set  of 
numerical  results,  including  simulation  results  for  a  directional 
slotted  Aloha  protocol  illustrating  the  need  for  inter-node 
coordination.  Finally,  we  provide  a  detailed  description  of  the 
directive  antenna  patterns  used  in  our  numerical  results. 

The  rest  of  this  paper  is  organized  as  follows.  Section 
II  outlines  related  work.  In  Section  III,  we  describe  our 
mm-wave  network  model.  In  Section  IV,  we  analyze  the 
statistics  of  interference  with  highly  directional  antennas, 
assuming  a  randomly  deployed  network  with  uncoordinated 
transmissions.  Section  V  presents  simulation  results  for  a 
directional  slotted  Aloha  protocol.  Finally,  conclusions  and 
future  work  are  discussed  in  Section  VI. 

II.  Related  Work 

There  is  a  significant  prior  (and  growing)  literature  on 
60  GHz  indoor  personal  area  network  design,  propagation 
modeling  and  link  measurements;  see  [11  ]— [  13]  and  the  refer¬ 
ences  therein.  While  MAC  design  ideas  developed  for  outdoor 
networks  may  be  applicable  in  this  domain  as  well,  additional 
considerations  such  as  blockage  and  the  coexistence  of  a 
number  of  different  applications  may  also  become  important 
in  the  indoor  setting. 

There  is  extensive  literature  on  directional  networking 
for  cellular,  broadband,  and  WiFi-based  multihop  wireless 
networks  operating  over  lower  frequency  bands  [14]— [19], 
Most  of  the  directional  networking  proposals  for  multihop 
wireless  networks  employ  a  combination  of  directional  and 
omnidirectional  communication  for  protocol  operation  and 
data  transfer.  These  approaches  exploit  the  broadcast  nature 
of  omnidirectional  transmissions  for  critical  control  message 
propagation,  thereby  avoiding  some  of  the  coordination  issues 
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that  arise  or  are  exacerbated  on  being  fully  directional  (such  as 
neighbor  discovery  and  deafness).  This  dual-mode  operation 
is  not  appropriate  for  the  mm-wave  mesh  networks,  where 
very  high  directionality  is  required  simply  to  achieve  a  reliable 
high  data  rate  link.  There  are  a  few  recent  proposals  on  fully 
directional  protocols  for  multihop  wireless  networks  [20]-[22]. 
However,  the  directivity  achievable  at  the  lower  frequency 
bands  considered  by  these  papers  is  much  smaller,  so  that 
the  focus  of  MAC  design  is  still  on  interference  management, 
unlike  the  pseudowired  abstraction  that  follows  from  the 
interference  analysis  in  our  paper. 

There  is  also  a  significant  body  of  previous  work  on 
modeling  the  statistics  of  interference  in  wireless  networks. 
Gupta  and  Kumar,  in  their  influential  contribution  on  wireless 
capacity  [9]  introduced  two  different  models  for  interference, 
i.e.,  the  “protocol  model”  and  “physical  model”.  The  protocol 
model  corresponds  to  the  collision  model  that  has  historically 
been  very  influential  in  MAC  design  for  packet-radio  networks 
[23],  while  the  physical  model  has  been  commonly  used 
to  model  co-channel  interference  [24]  in  cellular  networks. 
For  the  latter  application,  a  Gaussian  approximation  for  the 
aggregate  interference  [24],  [25],  motivated  by  the  Central 
Limit  Theorem  has  been  proposed.  However,  the  limits  of  the 
Gaussian  model  especially  for  channels  with  small  numbers  of 
users  and  highly  varying  channel  fading  gains  has  been  widely 
recognized  [26].  The  validity  of  the  simplifying  assumptions 
underlying  the  different  interference  models  was  explored  in 
[27],  which  also  provides  a  survey  of  the  models  used  in  MAC 
design.  The  authors  in  [27]  show  that  results  from  simple 
interference  models  similar  to  the  protocol  model  can  deviate 
significantly  from  a  more  accurate  model  based  on  the  SINR 
with  sum  interference  calculation  (i.e.,  the  physical  model). 

We  have  recently  learned  that  an  analytical  approach  similar 
to  ours  has  been  employed  in  the  information  theory  litera¬ 
ture  for  characterizing  outage  probabilities  in  omnidirectional 
networks;  see  [28]  and  the  references  therein.  The  approach 
is  conceptually  similar,  in  that  it  assumes  transmissions  are 
uncoordinated  and  distributed  according  to  a  Poisson  process 
on  a  plane,  and  that  the  interference  is  decomposed  into  “near” 
and  “far”  interferers.  Of  course,  the  details  of  the  analysis  for 
our  60  GHz  setting  are  quite  different  because  of  the  highly 
directional  nature  of  the  transmissions  and  the  exponential 
oxygen  absorption  loss. 

To  the  best  of  our  knowledge,  except  for  the  preliminary 
results  presented  in  the  conference  paper  [10]  that  we  expand 
upon  here,  the  present  paper  is  the  first  to  characterize  spatial 
interference  in  mm-wave  mesh  networks  and  to  quantitatively 
justify  a  pseudowired  abstraction  for  MAC  design.  We  have 
followed  up  on  these  MAC  design  guidelines  in  another 
conference  paper  [8],  where  we  present  a  MAC  protocol  that 
employs  memory  and  learning  to  address  deafness,  while 
exploiting  the  reduction  of  interference  between  simultaneous 
transmissions  to  take  a  reactive  rather  than  a  proactive 
approach  to  interference  management. 

III.  Network  Model 

We  consider  a  Poisson  distribution  of  nodes  over  a  large 
area  A  with  a  density  ps.  If  we  now  randomly  select  a  subset 
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Fig.  2.  Network  model  for  interference  analysis. 

of  nodes  as  transmitters,  the  distribution  of  transmitters  on  the 
area  of  interest  is  also  Poisson,  with  density  p  =  pspt  where  pt 
is  the  probability  that  the  selected  node  is  transmitting.  There¬ 
fore,  the  number  of  transmit  nodes  Nj-  over  a  deployment  area 
A  is  a  Poisson  random  variable  with  mean  pA. 

We  assume  that  each  node  can  communicate  with  at  most 
one  other  node  at  any  given  time  slot.  In  other  words,  we  do 
not  rely  on  advanced  physical  layer  capabilities  such  as  spatial 
multiplexing  or  multi-user  detection.  If  multiple  neighbors  are 
transmitting  to  the  same  receiver,  at  most  one  of  them  can  be 
successfully  decoded  by  the  receiver.  All  other  transmissions 
in  the  network  act  as  interference  for  the  receiver.  The  amount 
of  interference  depends  on  the  location  of  the  interferer  relative 
to  the  receiver,  and  the  radiation  patterns  of  the  antennas  at 
the  receiver  and  the  interferer.  We  assume  that  a  transmission 
is  successfully  decoded  by  the  receiver  if  the  total  signal  to 
interference  and  noise  ratio  (SINR)  is  above  a  given  threshold, 
say  (3.  Otherwise,  a  collision  occurs  and  the  transmission 
is  lost.  For  simplicity  of  analysis,  we  ignore  thermal  noise 
and  only  consider  interference  while  calculating  SINR. 
Similar  network  models  that  draw  upon  stochastic  geometry 
(e.g.,  randomizing  transmitter  locations  over  the  network 
deployment  area)  have  been  used  for  transmission  capacity  and 
interference  analysis  for  low  frequency  wireless  networks  [28]. 

A.  Nominal  link 

The  standard  Friis  transmission  equation  for  mm-waves 
gives  the  received  power  as  a  function  of  range  r  as 

PRx{r)  =  PTxGRxGTx(^)2e-ar  (1) 

where  I’tx  is  the  transmitted  power,  GRx,  Gtx  are  the  gains 
of  the  receive  and  transmit  antennas  respectively,  A  is  the 
wavelength,  and  a  is  the  attenuation  factor  due  to  absorption 
in  the  medium.  For  an  mm  wave  link  at  60  GHz,  A  =  5 
mm.  Parameter  a  captures  the  oxygen  absorption  loss  that 
can  be  as  high  as  16  dB/km  [29]  (for  oxygen  absorption 
of  16  dB/km,  a  =  0.0016  *  loge(10)  =  0.0037  per  meter). 
Since  lower  absorption  rate  leads  to  more  interference,  we 
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use  a  conservative  oxygen  absorption  value  of  10  dB/km 
(a  =  0.0023  per  meter)  in  our  numerical  results. 

We  assume  that  each  transmitter-receiver  pair  use 
beamsteering  to  maximize  the  signal  power  on  the  respective 
link  without  regard  to  the  interference.  Thus,  the  received 
power  at  a  reference  link  distance  of  Rq  (e.g.,  Rq  =100m  and 
200m  in  our  numerical  results)  is  denoted  by  Po  =  Prx(Ro), 
with  GTx  =  Grx  =  Gmax,  where  Gmax  denotes  the 
maximum  antenna  gain. 

Consider  the  link  budget  for  a  2  Gbps  line  of  sight 
(LoS)  link,  assuming  quadrature  phase  shift  keying  (QPSK) 
signaling,  a  desired  SINR  of  14  dB  (which  allows  for  uncoded 
QPSK  modulation  at  a  bit  error  rate  (BER)  of  10-9),  oxygen 
absorption  loss  of  10  dB/km,  and  10  mW  transmit  power.  For 
a  desired  link  range  of  Po  =  100m,  we  need  antenna  gains 
of  about  24  dBi  at  both  transmitter  and  receiver  in  order  to 
provide  a  10  dB  link  margin.  For  link  range  P0  =  200m, 
the  required  antenna  gains  at  the  transmitter  and  the  receiver 
increase  to  27  dBi.  These  link  parameters  are  used  as  a 
baseline  for  the  rest  of  the  paper. 

Our  setting  for  interference  analysis  is  as  follows.  Consider 
the  transmitter-receiver  pair  shown  in  Fig.  2.  Without  loss  of 
generality,  assume  that  the  receiver  is  located  at  the  origin 
and  is  communicating  with  the  transmitter  located  along  the 
X-axis,  at  a  distance  equal  to  the  reference  link  distance 
Po  while  undergoing  interference  from  other  concurrent 
transmissions.  The  other  Nt  —  1  interfering  transmitters  are 
transmitting  to  corresponding  receivers  located  at  randomly 
chosen  orientations. 


B.  Antenna  Pattern 


Directional  antennas  are  characterized  by  their  pattern  func¬ 
tions  that  measure  the  power  gain  G(0,  (f>)  over  the  spherical 
elevation  and  azimuthal  angle  coordinates  6  and  cj> .  The 
directivity  D  of  an  antenna  (which  equals  the  antenna  gain  for 
a  lossless  antenna)  is  the  ratio  of  the  maximum  power  density 
and  the  average  power  density  over  a  sphere.  It  is  given  by 


D  = 


An 

TP 


(2) 


where  f l  =  JJL0  JJ=_ w  Gn(9,  (j>)  sin  6  dO  d<f>  is  the  beam  area 
or  beam  solid  angle  of  the  antenna  in  steradians,  and  Gn(9,  <j>) 
is  the  normalized  power  pattern  of  the  antenna. 

Without  loss  of  generality,  we  assume  that  all  nodes  are  on 
the  same  horizontal  plane.  We  do  not  consider  variation  in 
beam  pattern  over  the  elevation  angle  6,  and  work  with  the 
normalized  2-dimensional  pattern 


5(0)  = 


G(  o,0) 


where  Gmax 


maxG(0,  <f>). 

<t> 


(3) 


We  define  the  azimuthal  beam  angle  of  the  antenna  as 

A0  =  f  5(0)#-  (4) 

J  (j)=—7T 


Fig.  3.  Approximating  a  circular  array  of  slot  antennas  as  a  uniform  linear 
array  of  flat-top  elements. 


antenna,  which  has  unit  gain  within  its  beam  angle  A (f>ft  and 
zero  gain  outside.  More  precisely. 


5(0)  = 


|0| 

otherwise. 


(5) 


In  3-dimensions,  the  flat-top  antenna  beam  is  assumed  to  be 
symmetric  about  the  beam-axis,  so  the  azimuthal  beam  angle 
and  the  elevation  beam  angle  are  equal. 

While  the  flat-top  antenna  is  a  useful  idealization,  practical 
directional  antenna  gains  have  a  more  complex  dependence  on 
the  azimuth  angle.  For  instance,  side  lobes  in  the  gain  function 
could  cause  significant  interference  even  in  directions  far  from 
the  antenna  boresight.  While  exact  computation  of  the  gain 
functions  of  practical  mm-wave  antenna  arrays  can  be  messy 
(because  array  elements  are  directional  themselves),  we  can 
obtain  useful  models  with  some  simplifying  assumptions.  In  a 
circular  array  of  directional  antenna  elements,  only  a  subset  of 
the  array  elements  with  significantly  overlapping  main  lobes 
contribute  to  the  antenna  gain  in  any  specific  direction.  The 
number  of  elements  in  this  subset  can  still  be  quite  large:  for 
instance,  if  the  diameter  of  the  overall  array  is  10  cm,  and  the 
inter-element  spacing  is  a  half-wavelength,  and  each  element 
has  an  azimuthal  beamwidth  of,  say,  120°,  we  can  have  on 
the  order  of  20  —  30  elements  with  overlapping  beams.  If  we 
assume  that  the  beams  of  these  elements  are  approximately 
flat  within  the  main  lobe,  and  neglect  the  curvature  of  their 
placement,  we  can  approximate  this  sub-array  as  a  uniform 
linear  array,  each  of  whose  elements  has  a  flat-top  response. 
In  particular,  we  obtain  the  following  normalized  power  gain 
pattern  for  an  N  element  linear  array  in  which  each  individual 
flat-top  element  has  beam  angle  (or  sector  size)  A <pft,  and  is 
placed  ^  apart: 

{1  sill3  (J-7T  sin  <t>)  1,1  ,  A <j>ft 

N 2  sin2(±7rsin  0)  ’  l  M-  2  >  (g) 

0,  otherwise. 


Fig.  3  illustrates  how  to  approximate  a  circular  array  with 
a  linear  array  of  flat-top  elements.  Fig.  4  shows  the  gain 
patterns  for  a  narrow  beam  flat-top  antenna  with  beam  angle 
14.4°  and  a  12-element  linear  array  of  flat-top  elements, 
each  of  sector  size  20°.  The  antenna  directivity  in  both  cases 
is  24  dBi.  The  corresponding  directivity  computations  are 
described  in  the  Appendix. 


IV.  Interference  Analysis 


One  idealization  that  proves  to  be  very  useful  for  our  inter-  We  now  analyze  the  probability  of  packet  loss  due  to 
ference  analysis  is  that  of  a  sectorized  “flat-top”  directional  interference.  We  consider  two  different  models  of  interference. 
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90°  90° 


(a)  Flat-top  antenna  (b)  Linear  array 


Fig.  4.  Gain  pattern  for  a  flat-top  antenna  and  a  linear  array  of  flat-top 
elements. 


In  the  protocol  model  of  interference,  a  packet  loss  occurs  if 
and  only  if  there  is  some  interfering  node  whose  signal  at  the 
receiver  exceeds  a  given  threshold.  Nodes  located  far  away 
from  the  receiver  cannot  cause  interference  in  this  model.  In 
the  physical  model,  a  packet  loss  occurs  when  the  total  inter¬ 
ference  from  all  nodes  exceeds  a  given  threshold.  Under  this 
model,  it  is  possible  for  the  sum  interference  from  a  number 
of  nodes  to  cause  packet  failure,  even  if  the  interfering  signal 
from  each  of  them  may  be  individually  too  weak  to  cause  a 
collision.  The  two  models  can  be  summarized  as  follows: 

Pr(collision)  = 

j Pr  ( max1<fe<^r_i  Pk  >  jPq),  (protocol  model) 

\Pr  (  Efc=Ti-1  pk  >  jPo),  (physical  model), 

where  Pk  is  the  power  at  the  receiver  of  the  signal  from  the  kth 
interferes  Since  pk  >  max1<i.<jvT_i  Pk,  it  follows 

that  Pr  (collision)  for  the  physical  model  is  lower-bounded 
by  the  corresponding  value  for  the  protocol  model. 

We  first  derive  an  expression  for  the  collision  probability  for 
a  “typical”  node  in  the  interior  of  a  large  network  under  the 
protocol  model  for  the  flat-top  antenna.  The  flat-top  antenna 
interference  analysis  under  the  protocol  model,  while  being 
simple  and  intuitive,  offers  valuable  insight  into  the  key  trends 
obtained  from  the  more  rigorous  models  presented  later  in  this 
section.  We  then  extend  the  analysis  to  general  directional 
antenna  models  that  account  for  side-lobes.  We  next  estimate 
Pr(collision)  with  the  physical  model. 

Our  analytical  approach  for  the  protocol  model  is  inspired 
by  the  analysis  of  localization  error  in  [30].  We  note  that 
our  analysis  applies  to  the  network  as  a  whole.  The  collision 
probability  for  a  typical  node  in  the  interior  of  a  large  network 
provides  a  pessimistic  estimate  of  the  average  collision  rate 
over  the  network,  since  it  ignores  edge  effects  (i.e.,  nodes  on 
the  edge  of  the  network  would  experience  less  interference 
than  those  in  the  interior). 

A.  Protocol  model  with  ideal  flat-top  antennas 

For  the  ideal  flat-top  antenna,  only  interferers  located  within 
the  boresight  of  the  receiver  can  cause  a  collision  (see  Fig.  5). 
Further,  a  transmitting  node  within  this  sector  causes  inter¬ 
ference  only  if  the  receiver  is  within  its  boresight,  which 
has  probability  q  =  since  the  potentially  interfering 

transmitter  is  sending  to  a  randomly  chosen  receiver. 


Fig.  5.  Protocol  model  with  flat-top  antennas. 


Fig.  6.  Interference  range  R^  for  Rq  =  100m,  from  (9). 


Let  Ri  be  the  interference  range,  i.e.,  the  maximum  distance 
an  interferer  can  be  from  the  receiver  and  still  cause  a  collision. 
Using  (1),  the  signal  and  interference  powers  are  evaluated  as 

PRx(R0)  =  P0  =  PTxGlax(^ye-aR°  (7) 

Pinter f(Ri)  =  Pt xGmax  (  /[n£>.  )  e 

=  P0§e-Q(fli-jRo)  (8) 

Ri 

where  we  used  Gtx  =  Grx  =  Grnax  for  the  antenna 
gains  assuming  that  the  interferer  and  the  receiver  are  within 
each  other’s  boresights.  We  set  Prx(Rq)  =  Po:  when  the 
transmitter  and  receiver  are  steered  towards  each  other,  this 
is  the  signal  power  designed  for  at  the  reference  distance  Rq. 
Using  the  collision  condition  for  obtaining  the  interference 
range  Ri,  Pinter f(Ri)  =  ^Po,  we  can  rewrite  (8)  as 


Ri 

if 


x(Ri-Ro)  _ 


=  A 


(9) 


which  determines  Ri  as  a  function  of  the  SINR  threshold  3. 
Fig.  6  plots  Ri  as  a  function  of  f3  for  an  mm-wave  mesh  link 
of  range  Rq  =  100m,  with  a  =  0.0023  per  meter. 

The  number  of  potentially  interfering  transmitters  is  there¬ 
fore  a  Poisson  random  variable  with  mean  p?  =  pA,,  where 
A,  =  (j>Rf.  The  probability  of  any  of  these  actually  causing 
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Fig.  7.  The  geometry  of  interference  with  directional  antennas. 


a  collision  is  q  =  2  ,  so  that  the  number  of  interferers  N, 
causing  a  collision  is  also  a  Poisson  random  variable,  with 
mean  g/r,  .  The  probability  of  a  collision  is  therefore  given  by 


Pr(collision)  =  Pr (iV*  >  0) 

(Aj.)2 

=  1  —  e 


=  1  -  e"w* 

_  ]_  _  g-P^c 


where  Ac 


(A0)2 

47T 


(3Rle-a{-Ri~Ro) . 


00) 

(11) 


From  (10),  Ac  =  cflnR2  can  be  interpreted  as  the  “expected 
interference  area”  around  the  receiver,  equal  to  the  region 
spanned  by  the  interference  range  ( tt  R2 )  scaled  down  by  q2, 
the  probability  that  the  receiver  and  the  interferer  are  within 
each  other’s  boresight. 

For  a  beam  angle  of  A cj>  =  ^  radians  (10°),  (3  =  15  dB, 
and  pi?Q  =  1  corresponding  to  roughly  ~pR2  ~  3  potential  in¬ 
terferers  within  communication  range  R0  =  100m,  and  50  po¬ 
tential  interferers  within  the  interference  range  Ri  =  400m  of 
each  receiver,  (10)  gives  an  estimate  of  Pr(collision)  ss  3.8%. 
The  substantially  low  collision  probability  in  this  example 
suggests  that  acceptable  MAC  performance  may  be  possible 
with  minimal  coordination  for  interference  management. 


B.  Protocol  model  with  general  directional  antennas 

We  now  generalize  (10)  for  a  general  directional  antenna 
pattern.  We  first  compute  the  probability  of  collision  due 
to  a  single  interferer  randomly  located  at  a  distance  R,  and 
angle  <j)\  relative  to  the  receiver  as  shown  in  Fig.  7.  The 
angle  02  represents  the  direction  of  the  interferer’s  beam 
relative  to  the  receiver.  We  model  <f>i,  02  as  independent 
and  uniformly  distributed  over  (— 7r,7r],  given  the  random 
orientation  of  the  interfering  transmitter  and  its  beam  relative 
to  the  desired  receiver.  Signal  power  is  still  given  by  (7),  and 
the  interference  power  is 

Pinterf(R)  =  (^) V“fl,  (12) 

where  we  used  GRx  =  Gmaxg{(j) i)  and  GTx  =  Gmaxg{^ 2). 
Using  (7),  we  can  rewrite  (12)  as 

Pinterf(R)  =  Po9(<t>l)g(<h)  ( jf)' .  (13) 


We  define  X  =  p,nlyXl<!  ^  the  normalized  interference  power 
due  to  a  transmitter  randomly  oriented  and  randomly  placed 
at  a  distance  R  from  the  receiver.  Therefore, 

X{fa,<h,R)  =  g{4>t)9{^)(^)2  e~a{R-Ro)  ■  (14) 

From  the  uniform  spatial  distribution  of  nodes,  the  conditional 
cumulative  distribution  function  (CDF)  FR(r)  of  the  random 
variable  R  representing  the  interferer  distance  (i.e.,  the 
probability  that  an  interfering  node  is  within  a  distance 
r  of  the  desired  receiver,  given  that  it  is  within  area  A) 
is  FR[r)  =  which  yields  the  conditional  probability 

density  function  (PDF)  fR(r)  as,  fR(r)  =  £(FR(r))  = 

The  distribution  of  X  is  determined  by  random  variables  R, 
(j> i,  and  02  that  are  mutually  independent.  Define  X,  as  the 
normalized  interference  caused  by  an  interferer  i:  X,  are 
independent  and  identically  distributed  (iid)  random  variables 
with  the  preceding  probability  distribution. 

Therefore,  the  probability  pr  that  an  interferer  would  cause 
a  collision  is 


Pc  =  P0X  >  i) 

=  II  j'  >  1) 

02——^ 


Ac 
A  ’ 


2n  r  1 

A  Air2 


dr  d(f>  i  d(f> 2 


(15) 


where  l(-)  is  the  indicator  function  that  takes  the  value  1  when 
its  argument  is  true,  and  0  otherwise,  and 

=  Jl  f  i(g(ti)g(fo)^e-^-R^  >  i) 

(r,0iGA)  02=—™ 

—rdr  d(j)\  dtp 2.  (16) 

2n 


We  now  consider  Nt  —  1  interferers  placed  randomly  in  the 
area  A.  Each  interferer  has  a  collision  probability  pc  with  the 
receiver  given  by  (15).  A  collision  occurs  if  at  least  one  of 
these  interferers  causes  a  collision.  The  probability  of  collision 
is  given  by 


Pr(collision) 


1  -  (1  -  Pc)^-1 

Ac 

~A 


1  _  lim  (1  _ 

00  V  A  ’ 


1  -e~pAc. 


(17) 

(18) 


Note  that  (18)  has  an  identical  form  to  (10),  and  that  the 
collision  probability  depends  on  the  antenna  pattern  only 
through  Ac.  Thus,  for  the  protocol  model,  we  can  restrict 
attention  to  an  equivalent  flat-top  model  in  the  azimuthal 
plane  where  the  general  directional  antenna  of  each  node  is 
replaced  by  a  flat-top  antenna  whose  beam  angle  A (f>eq  can 
be  calculated  by  evaluating  Ac  using  (16),  and  plugging  into 
(11)  to  obtain  A (j)eq  =  'jppr- ^Ri~ R°\  For  instance,  a 
6-element  linear  array  of  flat-top  antennas  of  sector  size  120° 
and  half-wavelength  spacing  has  an  equivalent  flat-top  beam 
angle  of  about  23°  for  a  =  10  dB/km,  Rq  =  100m,  and 
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Fig.  8.  Equivalent  flat-top  beam  Fig.  9.  Equivalent  flat-top  beam 
angle  (link  range  =  100m).  angle  (link  range  =  200m). 
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Fig.  10.  Interference  with  the  physical  model:  a  hybrid  approach. 


(3  =  15  dB.  The  algebraic  azimuthal  beam  angle  for  the  same 
linear  array,  given  by  (4)  is  about  20°.  Figs.  8  and  9  show  the 
equivalent  flat-top  beam  angles  for  linear  arrays  of  different 
numbers  of  flat-top  elements  of  sector  size  120°  and  half¬ 
wavelength  spacing,  for  communication  ranges  of  100m  and 
200m.  As  seen  from  the  figure,  the  equivalent  flat-top  beam 
angle  is  numerically  close  to  the  algebraic  beam  angle  given 
by  (4)  and  does  not  vary  much  with  the  SINR  threshold  (3. 


C.  Physical  Model 

In  the  protocol  model,  only  nodes  located  within  a  bounded 
distance  from  the  receiver  are  capable  of  causing  a  collision. 
On  the  other  hand,  under  the  physical  model,  interfering 
signals  from  a  large  number  of  far-away  transmitters  could, 
in  principle,  sum  up  at  the  receiver  to  cause  packet  reception 
failure.  It  turns  out  that  the  statistical  properties  of  interference 
random  variable  X(<f>x,  fa,  R)  (defined  in  (14))  are  such  that 
the  tail  probabilities  do  not  die  fast  enough,  leading  to  a 
probability  distribution  skewed  away  from  the  small  mean 
value  in  the  positive  direction,  with  a  high  variance.  We  find 
that  the  traditional  upper  bound  techniques  such  as  Chemoff 
and  Markov-type  bounds  do  not  work  well  when  applied 
directly  to  characterize  the  sum  interference  over  a  large  area. 


Fig.  11.  Pr  ( Xfar  >  jjAX)  for  Ro  =  100m.  from  (21). 


Furthermore,  applying  the  Central  Limit  Theorem  (CLT)  on 
the  sum  interference  does  not  yield  reasonable  Pr(collision) 
estimates  for  typical  node  densities  of  interest. 

While  the  protocol  model  collision  probability  serves  as  a 
lower  bound  for  the  collision  probability  under  the  physical 
model,  we  separately  consider  the  effects  of  “far-away” 
and  “nearby”  interferers  (say,  within  a  circular  region  C Rlh 
of  radius  Rth  centered  at  the  receiver)  in  order  to  obtain 
an  approximation  and  an  upper  bound  for  the  collision 
probability.  We  bound  the  effect  of  far-away  interferers 
by  applying  the  Markov  inequality.  We  employ  techniques 
similar  to  the  protocol  model  to  capture  the  effect  of  dominant 
nearby  interferers  (i.e.,  that  can  individually  cause  a  packet 
loss  at  the  receiver),  while  the  aggregate  interference  from 
the  non-dominant  nearby  interferers  is  characterized  via  the 
CLT  to  obtain  an  approximation,  and  the  Chernoff  bound 
to  compute  an  upper  bound  for  the  collision  probability.  We 
find  that  the  CLT-based  approximation  provides  an  accurate 
estimate  for  the  collision  probability  over  a  wide  range  of 
parameter  choices  even  when  the  bounds  are  not  necessarily 
very  tight.  Fig.  10  illustrates  the  setting. 

Let  r,  be  the  distance  of  the  ith  interferer  from  the  receiver. 
We  write  the  total  normalized  interference  power  as  the  sum 
of  two  contributions  Xnear  and  Xfar,  defined  as: 

Xnear  =  XU  and  Xfar  4  ^  X,;  (19) 

{i:ri<Rth}  {i:r  i  >Rth  } 

where  Rtu  is  a  suitable  large  distance.  Then  we  have 
Pr  (collision)  =  Pr  (Xnear  +  Xfar> 

<  Pr  (Xfar  >  AA^)  +  Pr  ( Xnear  >  Xt), 

(20) 

where  A  A  +  A t  =  4.  Xt  and  A  A  are  defined  as  thresh¬ 
olds  for  the  nearby  and  far  away  interferers.  If  we  choose 
AA'  -C  A t,  the  second  probability  Pr  ( Xnear  >  Xt)  ~ 
Pr  (Xnear  >  jj).  Our  approach  in  (20)  is  to  bound  away 
the  effect  of  far-away  interferers  and  to  compute  analytical 
estimates  for  the  effect  of  nearby  interferers,  which  are  ba¬ 
sically  the  cause  of  packet  loss.  However,  in  order  to  do 
this  effectively  and  obtain  close  Pr(collision)  estimates,  we 
must  choose  radius  Rth  large  enough  such  that  the  far-away 


interferers  are  insignificant  enough.  We  choose  Rth  =  50f?o 
and  AXdB  =  —fidB  —  30  dB,  which  ensures  that  Pr  ( Xfar  > 
AX)  <  0.02  Pr  (Xnear  >  XT)  for  all  SINR  values  /?  in  our 
setting. 

We  can  bound  the  first  term  in  (20)  using  the  Markov 
inequality: 

Pv(Xfar  >AX)  <  (21) 


The  expectation  E  [Xfar\  in  (21)  can  be  readily  evaluated  as: 


E^]  -  £ 


I  f  g{<h)9{<h)^e 


0  -oc(r-R0) 


r=Rth 


rdr  d<f>i  d(f> 2 


(A 

2tt 


pR20eaRo 


-dr, 


(22) 


r=Rth 


where  Ap  is  the  antenna’s  azimuthal  beam  angle  as  defined  in 
(4).  Fig.  11  plots  Pr  ( Xfar  >  AX)  as  a  function  of  (3,  using 
(21),  for  R0  =  100m  assuming  flat-top  antenna  model  and 
different  interferer  node  densities.  We  will  see  in  Section  IV-D 
that  these  values  are  more  than  an  order  of  magnitude  lower 
than  the  estimated  collision  probabilities  for  these  settings. 
The  effect  of  nearby  interferers:  We  now  estimate 
Pr  ( Xnear  >  Xt)  in  (20),  where  we  first  present  a  simple 
approach  for  the  flat-top  antenna  model  and  then  consider 
general  directional  antennas. 

Each  transmitter  i  inside  the  circle  C Rth  causes  a  random 
amount  of  interference  Xi[<pi,  fa,  R)-  The  interference 
random  variables  {2Q}  are  independent  and  identically 
distributed  (iid),  and  are  also  independent  of  the  number  of 
interferers  inside  Cnth,  say  Nth  ( Nth  is  Poisson  distributed 
with  mean  and  variance  of  7 th  =  P^R2^,  since  conditioned 
on  Nth,  the  location  of  any  given  interferer  is  iid  and  uniform 
inside  C a jh .  Note  the  two  conflicting  effects  of  distance 
R:  the  probability  of  there  being  an  interferer  increases 
with  distance  R  from  the  receiver  (/^j(r)  =  2N  ),  but  the 
resulting  interference  (given  by  (14))  is  drastically  low  for 
each  interferer  located  at  larger  R. 

Flat-top  antenna  model:  For  a  flat-top  antenna,  the  interfer¬ 
ence  threshold  Xt  is  exceeded  by  the  nearby  interferers  under 
two  events,  A  and  B,  defined  as  follows: 

A:  there  is  at  least  one  interferer  that  individually  causes 
interference  >  Xt-  Note  that  such  an  interferer  must  be  within 
a  bounded  distance,  say  lip  of  the  receiver  (lip  ss  II,  defined 
in  (9)  for  AX  «  XT-)- 

B:  the  interference  beyond  Up  sums  up  to  >  Xt- 

Events  A  and  B  are  independent  events  as  they  correspond 
to  transmitters  located  in  two  disjoint  spatial  regions  (illus¬ 
trated  in  Fig.  12),  the  statistics  of  which  are  modeled  by 
two  independent  Poisson  processes:  A  results  from  interferers 
located  within  range  Rp  (circle  Cp )  around  the  receiver, 
whereas  B  depends  on  the  interferers  located  within  the 
annular  region  between  Rp  and  Rth-  Therefore, 


Pr  [Xnear  >  XT)  =  Pr(*4)  +  Pr(B)  -  Pr(*4)  Pr (B).  (23) 
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Fig.  12.  Flat-top  antenna:  estimation  of  Pr  ( Xnear  >  Xt)- 


Consider  event  A:  let  random  variable  Y(<f>  1,  (In  -  H  \ )  denote 
the  interference  caused  by  an  interferer  located  within  circle 
Cp.  Here,  fR1[r)  =  for  r  <  Rp,  and  <j>i  and  <f>2  are  uni¬ 
formly  distributed  over  ( — 7r,  7t]  as  earlier.  Note  that  in  terms 
of  the  general  interference  random  variable  X((j>\,  fa,  R), 
Y  =  X{<l>1,<i>2,R\R<RA). 

The  analysis  for  Pri^tj  is  analogous  to  that  in  Section  IV-B 
for  the  protocol  model,  where  we  use  (15)  with  A  =  "  I l‘p  to 
obtain  the  probability  of  collision  pc  for  an  interferer  placed 
at  random  within  Cra.  Conditioned  on  Np,  the  number  of 
interferers  inside  Cra,  we  can  estimate  the  probability  that 
at  least  one  interferer  among  all  interferers  in  Cra  causes 
interference  >  XT  as 


Pr  ( max  Pi  >  Xt  \  Np  =  fc)  =  1  —  (l  —  pc)k-  (24) 


Pr  ( max,  V)  >  Xt)  can  then  be  computed  by  averaging  out 
Np,  using  the  fact  that  Np  is  Poisson  distributed  with  mean 
7*4  =  pTrRp.  We  can  also  substitute  7 p  in  place  of  k  in  (24) 
to  obtain  a  close  approximation. 

For  event  B,  let  Z(4>\ ,  (in,  II2)  denote  the  interference 
random  variable  corresponding  to  interferers  that  lie  within 
the  annular  region  between  Rp  and  Rth-  We  have  (r)  = 
for  Ra  <  r  -  Rth ■  Note  that  2  = 
X[(f>i,(j)2,R  |  Ra  <  R  <  Rth)-  We  apply  the  CFT  to 
estimate  Pr(JT  Z,  >  Xt)-  Conditioned  on  Arg,  the  number 
of  interferers  within  the  annular  region  between  Rp  and  Rth, 
we  have 

rr(S^^^I^  =  ‘)-Q(7=g!).  PS 


where  Q(  )  is  the  complementary  cumulative  distribution 
function  for  a  Normal  distribution,  given  by  Q(x)  = 
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f ^  e  2  dx.  The  mean  E [Z\  is  given  by 


E[Z\  = 


rRth 


Jf  g(Mg(^e 


0  -a(r-R0) 


r=RA  0i, 02=-ir 


fR2(r)^dr  dcj) i  dcf>2 


(A^y 


Rl 


2n  {R2th  -  R\) 


oaR0 


fRth 


dr.  (26) 


r=Rj i 


Var (Z)  can  be  evaluated  by  first  computing  E \Z2]  and  then 
using  Var(Z)  =  E [Z2]  —  (E [Z])~ . 

The  required  probability  estimate  Pr  ( Xnear  >  Xt)  can 
be  computed  by  averaging  out  (25)  over  NB.  Given  the  large 
expected  number  of  nodes  7 B  in  the  annular  region  between 
Ra  and  Rth,  given  by  =  pn(R2th  -  R2A),  (e.g.,  for  p  =  1 
and  Rth  =  50f?o,  78  =  P^iRth  ~  Ra)  ~  7-8  •  103R2)  we  can 
apply  the  CLT  with  mean  and  variance  computed  by  replacing 
Nb  with  E(Nb)  =  7b,  so  that  we  obtain 


Pr  (  ^2  Zi  >  Xt)  =  Q  ( 

i 


XT  -  7bE[^]  \ 
\J 7eVar(Z) 


(27) 


In  order  to  obtain  an  upper  bound  on  Pr  ( ZZ,  %i  X  Xt), 
we  apply  the  Chernoff  bound  as  follows: 

Pr  {Y,Zi>  Xt)  <  (E[eSxT]j  e~a,  (28) 

i 

where  the  Chernoff  bound  is  derived  by  minimizing  this  family 

/  s_g±_  \  "re  _ 

of  bounds  over  s  >  0.  Since  mins  >0  E[e  A>]  e  s  has  to 
be  obtained  via  an  exhaustive  search  over  a  sufficiently  large 
range  of  values  of  s  >  0,  the  calculation  of  the  Chernoff  bound 
turns  out  to  be  computationally  expensive. 

The  flat-top  antenna  analysis  presented  above  does  not  apply 
to  general  directional  antennas  because  the  independence  of 
events  A  and  B  does  not  hold:  for  general  directional  antennas, 
even  nodes  within  range  RA  of  the  receiver  can  cause  non-zero 
interference  less  than  Xt  (depending  on  the  relative  antenna 
orientations  and  gain  patterns)  due  to  which,  event  B  can 
include  transmitters  within  C  nA .  This  is  unlike  the  flat-top 
antenna  case  where  the  antenna  gains  at  both  the  interferer 
and  the  receiver  are  either  1  or  0,  leading  to  either  interference 
>  Xt  if  both  nodes  are  within  each  other’s  antenna  boresights, 
or  0  otherwise  (which  leads  to  the  independence  of  events  A 
and  B). 

General  directional  antennas:  We  now  extend  the  preceding 
analysis  to  general  directional  antenna  patterns,  where  we 
derive  a  CLT-based  approximation  and  obtain  an  upper  bound 
for  Pr  (X near  >  Xt)  ■  We  use  the  following  expression  to 
estimate  Pr  ( Xnear  >  Xt)' 

Pr  (Xnear  >  XT)  =  Pr  (maxi,  >  XT)  + 

(1  -  Pr  ( maxi,  >  XT))  Pr  ( ^  >  XT  \  X,  <  XT  V  i) . 


The  first  term  in  the  right-hand  side  (RHS)  of  (29)  captures  the 
probability  of  cases  when  the  interference  from  at  least  one 
interferer  exceeds  Xt  at  the  intended  receiver.  The  second 


term  captures  the  cases  when  there  is  no  single  dominant 
interferer  i  with  X)  >  Xt,  but  the  sum  interference  from 
all  the  interferers  exceeds  Xt-  Unlike  Xt,  the  conditional 
random  variable  (Xi  \  X,  <  Xt)  is  bounded,  its  variance  is 
much  smaller  than  that  of  Xi:  and  the  CLT  is  better  behaved 
for  settings  with  moderate  node  densities.  Therefore,  as  we 
describe  in  this  section,  (29)  is  amenable  to  obtaining  a 
close  approximation  via  the  CLT  and  an  upper-bound  via  the 
Chernoff  bound. 

We  evaluate  the  first  term  Pr  ( max,  X.t  >  Xt)  of  (29) 
via  the  protocol  model  approach  (Section  IV-B),  in  a  manner 
analogous  to  that  described  for  obtaining  Pr(_4)  in  the  flat-top 
antenna  model  case.  We  estimate  Pr  ( JT  X,  >  Xt  \  Xi  < 
Xt  V  i)  in  the  second  term  in  the  RHS  of  (29)  by  applying 
the  CLT  to  finally  obtain: 

Pr  ( Y,  Xi  >  XT  I  Xi  <XTVi)  = 

i 

'XT  -  7t/,  E[A"i  |  X, 


Q(: 


<  Xt\  \ 


where  E  [Xt  \ 
E[Xi  |  X, 

rRth 


^7thVar(Xi  |  X,  <  XT) 
Xi  <  Xt]  is  given  by 
;  <  Xt]  = 
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dr  d<pi  dcj) 2,  (31) 


where 


Pr  (Xi  <  XT)  = 

fth  ff  i(g(ch)g(<t>2)^e-a(r-Ro)  <  Xt) 
r—0  01,02  =  —  7T 

fR(r)fdr  d(f>  1  # 2.  (32) 

l(-)  above  is  the  indicator  function.  For  calculating 
Var(X,;  |  XL  <  Xt),  we  first  obtain  E [X?  |  X,;  <  Xt]  using 
a  similar  approach  as  shown  for  E [X,  |  X,  <  Xt],  and  use 
Var(X,;  |  X,  <  XT)  =  E[X?  |  X,  <  XT]  -  (E[X4  |  X,  < 
XT])2. 

We  show  in  Section  IV-D  that  (29)  evaluated  using  the  CLT, 
when  plugged  into  (20)  provides  excellent  estimates  of  the 
interference  probability  for  all  scenarios  of  interest  in  terms 
of  node  densities  as  well  the  antenna  patterns. 

In  order  to  obtain  an  upper  bound  on  Pr  (Xnear  >  Xt), 
we  first  calculate  the  Chernoff  bound  for  Pr  ( ZZ,  Xi  > 
XT  I  Xi  <  XT  V  i)  as: 

Pr  (^X,  >  XT  I  X,  <  XT  V  i)  < 

i 

min  (E[e  |  X*  <  XT]J  e  s.  (33) 

We  then  plug  this  bound  in  (29)  to  obtain  the  desired  upper 
bound.  We  present  the  Chernoff  bound  results  for  all  the 
examples  considered  in  the  next  section. 
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D.  Interference  statistics 

We  now  employ  our  analytical  models  to  estimate  collision 
probabilities  for  mm-wave  mesh  networks  considering  design 
parameters  described  in  Section  III.  In  particular,  we  study  the 
dependence  of  interference  statistics  on  antenna  directivities 
obtained  from  the  link  budget  analysis  for  link  ranges  of  100m 
and  200m,  and  spatial  deployment  density  of  mesh  nodes. 

The  collision  probabilities  for  the  protocol  model  are  com¬ 
puted  from  (18)  and  serve  as  lower-bounds  for  the  correspond¬ 
ing  collision  probabilities  from  the  physical  model.  We  use 
(20)  for  the  physical  model.  The  first  term  in  (20)  is  computed 
from  (21)  with  Rth  =  50Rq  and  A X^b  =  —/3dB  —  30  dB.  The 
second  term  in  (20)  is  estimated  via  the  following  methods: 
(a)  using  (29)  (we  also  include  results  from  (23)  for  the  flat- 
top  antenna  model  cases  for  comparison),  and  (b)  the  use  of 
Monte-Carlo  simulation,  which  also  serves  as  a  baseline.  In 
addition,  by  applying  the  Chernoff  bound  in  place  of  the  CLT 
for  calculating  the  second  term  in  (20),  we  obtain  an  upper 
bound  on  Pr(collision)  for  all  our  example  settings. 

We  first  consider  a  set  of  randomly  located/oriented  60  GHz 
outdoor  point-to-point  links.  The  transceiver  hardware  for  such 
line-of-sight  links  is  commercially  available  (e.g.,  [31]).  In  or¬ 
der  to  model  this  setting,  we  work  with  moderate  spatial  node 
densities  that  do  not  necessarily  ensure  a  connected  network. 

Figs.  13(a)  and  13(b)  show  the  collision  probabilities 
corresponding  to  different  values  of  the  SINR  target  j3,  for 
link  range  Rq  =  100m,  with  ideal  flat-top  antenna  and  with 
a  linear  array  of  flat-top  elements,  respectively,  both  with  a 
directivity  of  24  dBi  (required  to  close  the  100m  link;  see 
Section  III-A).  The  beam  angle  for  a  flat-top  antenna  with 
this  directivity  is  14.4°,  which  can  be  obtained  from  (2) 
with  the  flat-top  antenna  beam  pattern  (see  Appendix).  The 
linear  array  comprises  12  flat-top  elements,  each  of  sector 
width  20°,  placed  half  wavelength  apart.  The  azimuthal 
beam  angle  for  the  array  (given  by  (4))  is  8.7°,  whereas 
the  azimuthal  and  vertical  half-power  beamwidths  are  8.5° 
and  20°,  respectively.  The  spatial  node  density  is  such  that 
ttpRq  =  7 r,  and  oxygen  absorption  is  10  dB/km. 

We  observe  that  for  a  flat-top  antenna,  when  the  desired 
SINR  (3  increases  beyond  about  15  dB,  the  probability  of 
collision  approaches  10%,  whereas  for  the  linear  array,  owing 
to  the  lower  azimuthal  beam  angle,  the  probability  of  collision 
is  relatively  low,  and  approaches  10%  for  /3  =  21  dB. 
To  put  these  results  in  perspective,  note  that  even  a  single 
interferer  within  the  interference  range  of  the  receiver  would 
be  enough  to  cause  a  collision  in  an  omnidirectional  setting. 
We  observe  that  the  collision  probability  estimates  under  the 
physical  model  using  the  CLT  approach,  and  those  obtained 
with  Monte-Carlo  simulations  of  the  nearby  interferers  are 
very  close  to  each  other  and  to  the  corresponding  protocol 
model  estimates.  Moreover,  for  the  flat-top  antenna  model, 
the  physical  model  Pr(collision)  estimates  in  (20)  using  (23) 
to  account  for  the  effect  of  nearby  interferers  are  very  close  to 
those  obtained  via  (29).  We  also  observe  that  the  Pr(collision) 
upper-bounds  obtained  via  the  Chernoff  bound  are  loose  for 
this  setting,  especially  at  high  SINR  values. 

Figs.  14(a)  and  14(b)  show  the  interference  loss  probabilities 
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Fig.  13.  Collision  loss  probability  for  Rq  =  100m,  npR^  =  7r. 
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Fig.  14.  Collision  loss  probability  for  Rq  =  200m,  npR^  =  7r. 


for  ideal  flat-top  antenna  and  linear  array  of  flat-top  elements 
for  Rq  =  200m.  The  longer  link  range  of  200m  requires 
antenna  directivity  of  about  27  dBi,  which  corresponds  to  a 
flat-top  antenna  of  beam  angle  10.2°,  and  a  linear  array  of 
26  flat-top  elements  spaced  half  wavelength  apart,  each  with 
sector  width  of  20°.  The  azimuthal  half-power  beamwidth  of 
the  linear  array  is  about  4°.  The  other  parameters  are  the  same 
as  the  preceding  case  of  Rq  =  100m.  From  the  plots,  it  is 
clear  that  for  200m  link  range,  the  higher  antenna  directivity 
reduces  the  interference  loss  probability  even  further  relative 
to  that  for  100m  links.  Also,  in  this  highly  directional  setting, 
Pr(collision)  estimates  obtained  via  all  the  models  are  very 
close  to  each  other,  and  Pr(collision)  upper  bounds  obtained 
via  the  Chernoff  bound  are  relatively  tight.  Fig.  15  illustrates 
the  reduction  in  the  collision  probability  with  increasing  link 
range,  where  we  plot  Pr(collision)  for  flat-top  antennas  with 
different  beam  angles  corresponding  to  link  ranges  varying 


Fig.  15.  Collision  loss  probability  for  the  flat-top  antenna  model  with 
increasing  link  range  Rq  and  n pR^  =  n. 
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Fig.  16.  Collision  loss  probability  for  Rq  =  100m,  TtpRf  =  5.2. 


Fig.  17.  Collision  loss  probability  for  Rq  =  200m,  npR^  =  4.4. 


from  100m  to  350m,  with  the  other  link  budget  parameters 
same  as  in  Section  III-A. 

The  key  observation  from  the  above  results  is  that  the 
packet  loss  due  to  interference  is  expected  to  be  very  low  for 
mm-wave  mesh  networks  with  the  node  density  and  SINR 
parameters  of  interest  to  us.  Furthermore,  the  probability 
of  collision  under  the  physical  model  does  not  appreciably 
differ  from  the  protocol  model,  particularly  for  collision 
probabilities  of  10%  or  less,  implying  that  the  protocol  model 
estimates  work  well  in  this  desired  regime  of  operation. 

We  now  consider  higher  mesh  node  deployment  densities 
that  guarantee  a  connected  network  with  a  high  probability. 
Note  that  the  capability  to  sense  and  beam-steer  towards 
a  neighbor  in  any  azimuthal  direction  implies  that  the 
connectivity  versus  spatial  node  density  results  for 
omnidirectional  networks  apply  in  this  setting.  In  order 
to  calculate  the  required  node  densities,  we  employ  the 
analytical  model  and  connectivity  results  presented  in  [32]. 
Ignoring  border  effects  (which  is  reasonable  for  large 
networks)  we  obtain  that  for  Rq  =  100m,  npR^  =  5.2 
ensures  a  connected  network  with  a  high  probability  of  0.99, 
whereas  for  Rq  =  200m,  the  required  value  is  irpR g  =  4.4. 
Note  that  p  is  the  spatial  density  of  transmit  nodes,  which  is 
assumed  to  be  less  than  half  of  that  for  all  nodes.  Figs.  16(a) 
and  16(b)  show  the  collision  probabilities  for  the  ideal  flat-top 
antenna  and  for  the  linear  array  respectively,  for  link  range 
Ro  =  100m.  The  Pr(collision)  estimates  for  the  physical 
model  are  greater  than  those  obtained  under  the  protocol 
model  for  higher  values  of  the  desired  SINR  (e.g.,  >  18dB). 
The  reason  for  this  difference  is  that  the  increasing  effect  of 
the  sum  interference  (with  higher  node  densities)  that  leads 
to  more  collisions  for  higher  SINR  is  not  captured  by  the 
protocol  model.  We  also  observe  that  the  Pr(collision)  upper 
bounds  obtained  via  the  Chernoff  bound  are  very  loose  for  this 
case,  especially  for  higher  SINR  targets.  Figs.  17(a)  and  17(b) 
plot  the  interference  loss  probabilities  for  ideal  flat-top  antenna 
and  linear  array  of  flat-top  elements  for  Rq  =  200m.  Although 
the  probability  of  collision  in  each  case  is  comparatively 
higher  than  the  corresponding  setting  with  a  lower  node 
density  (Figs.  14(a)  and  14(b)),  the  previous  observation  of 
substantially  low  interference  loss  probabilities  relative  to 
omnidirectional  communication  settings  still  applies. 


E.  A  pseudowired  abstraction 

The  preceding  interference  analysis  results  show  that  for 
mm-wave  mesh  networks  with  moderate  to  high  node  den¬ 
sities,  even  with  uncoordinated  transmissions,  the  collision 
probabilities  are  very  small  for  the  parameters  corresponding 
to  the  reference  link  budget  described  in  Section  III-A.  This 
observation  points  to  a  very  different  paradigm  from  that 
normally  encountered  in  wireless  networks,  where  interference 
substantially  limits  spatial  reuse.  Furthermore,  the  probability 
of  collision  under  the  protocol  model  is  very  close  to  that  with 
the  physical  model  in  the  desired  regime  of  operation  with 
lower  collision  probabilities  (e.g.  <  10%)  and  for  the  purpose 
of  interference  analysis,  details  of  antenna  beam-patterns  can 
be  abstracted  away  using  the  notion  of  equivalent  flat-top 
beam  angles.  This  suggests  that  MAC  protocols  for  mm-wave 
mesh  networks  need  not  be  focused  on  minimizing  packet  loss 
from  interference.  Instead,  the  focus  of  MAC  design  should  be 
on  overcoming  the  challenging  problem  of  scheduling  under 
deafness  constraints.  We  conclude  that  the  MAC  designer  can 
use  the  following  pseudowired  abstraction  as  a  starting  point: 

1)  No  interference.  Transmissions  between  two  distinct 
pairs  of  nodes  are  unlikely  to  interfere  with  each  other, 
and  interference  can  be  largely  ignored  in  MAC  design. 

2)  Half-duplex  constraint.  Each  node  can  either  transmit  or 
receive  at  any  given  time,  but  not  both. 

V.  MAC  Simulations 

The  results  in  Section  IV  indicate  that  we  can  essentially 
ignore  interference  and  the  details  of  antenna  patterns  in  the 
design  of  medium  access  control  for  mm-wave  mesh  networks, 
although  their  effect  on  performance  must  be  quantified.  This 
is  a  distinct  advantage  over  directional  networking  at  lower 
carrier  frequencies,  where  it  is  difficult  to  achieve  such  high 
link  directivities.  On  the  other  hand,  the  high  directivity  of 
mm-wave  links  also  implies  that  mm-wave  mesh  nodes  are 
essentially  deaf  to  their  neighbors,  which  poses  a  challenge  in 
coordinating  transmissions. 

To  verify  the  pseudowired  abstraction  and  highlight  its  ef¬ 
fects  on  network  performance,  we  simulate  a  naive  directional 
slotted  Aloha  protocol.  Prior  studies  on  slotted  Aloha  with 
directional  communication  include  [20],  [33],  [34],  but  unlike 
these  papers,  our  goal  is  to  examine  the  relative  effects  of 
interference  and  deafness  on  performance.  We  note  that  far 
better  performance  can  be  obtained  using  more  sophisticated 
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MAC  designs.  Our  own  related  publication  [8]  proposes  one 
such  MAC  protocol. 

Our  approach  to  compare  the  packet  loss  because  of  inter¬ 
ference  and  lack  of  transmit-receive  coordination  for  slotted 
Aloha  is  the  following.  We  obtain  the  number  of  packets  that 
are  not  decoded  correctly  by  the  intended  receivers  because 
of  external  interference  (which  depends  on  the  corresponding 
SINR  at  each  intended  receiver).  We  define  transmit-receive 
coordination  loss  as  the  packet  loss  either  because  the  intended 
receiver  was  beamformed  in  a  direction  other  than  that  towards 
the  transmit  node  under  consideration  (leading  to  poor  antenna 
gain  in  the  direction  of  the  transmit  node),  or  because  the 
intended  receiver  was  also  transmitting  a  packet  during  the 
slot.  We  compare  the  ratio  of  these  losses  with  respect  to  the 
total  number  of  packet  transmissions. 

Simulation  set-up:  We  consider  random  network  topologies 
with  25  or  50  nodes  spread  over  a  500m  x  500m  flat  terrain. 
Every  node  initiates  one  constant  bit  rate  (CBR)  flow  to  each 
of  its  neighbors,  at  a  large  data  rate  (1  Gbps)  enough  to 
saturate  the  link.  Under  the  directional  slotted  Aloha  protocol, 
whenever  a  node  has  a  new  packet  to  transmit,  it  beamforms 
towards  the  direction  of  the  intended  receiver  and  transmits 
the  packet  in  the  next  slot.  If  the  node  does  not  receive  an 
ACK,  it  attempts  to  retransmit  the  packet  with  a  probability 
Pretx  over  the  next  slots.  A  node  returns  to  the  unbacklogged 
state  after  every  successful  packet  transmission. 

We  consider  a  sectorized  antenna  design  with  each  sector 
covered  by  an  array  of  high-gain  antenna  elements.  For  con¬ 
creteness,  we  consider  endfire  patch  antenna  elements  that  can 
be  realized  as  patterns  of  metal  on  a  circuit  board  [35].  Each 
element  has  azimuthal  and  vertical  halfpower  beamwidths  of 
about  33°  and  45°,  respectively,  which  leads  to  a  directivity 
of  14  dBi.  The  element  gain  patterns  are  shown  in  Fig.  18. 
In  order  to  obtain  the  desired  array  directivity  of  24  dBi,  we 
need  a  linear  array  of  25  elements  spaced  half  wavelength 
apart.  The  azimuthal  half-power  beamwidth  of  the  resulting 
array  is  about  4°. 

We  use  the  QualNet  network  simulator  [36]  for  our 
simulations.  We  have  modified  the  QualNet  PHY  and 
Antenna  modules  to  model  propagation  in  the  mm-wave  band 
and  our  link  budget  design. 

Results:  Figs.  19(a)  and  19(b)  show  the  interference  loss 
and  the  transmit-receive  coordination  failure  loss  for  25  and 
50  node  random  topologies  as  a  function  of  slotted  Aloha 


retransmit  probability  pretx ■  The  loss  figures  are  relative 
to  the  total  number  of  packet  transmissions.  Fig.  19(c) 
shows  the  corresponding  aggregate  network  throughput. 
Comparing  Figs.  19(a)  and  19(b),  we  observe  that  the 
interference  loss  is  about  an  order  of  magnitude  lower 
than  the  coordination  loss  in  this  setting,  which  indicates 
that  transmit-receive  coordination  is  a  significantly  bigger 
challenge  than  interference.  As  the  retransmission  probability 
increases,  nodes  become  more  aggressive  in  attempting 
transmissions  after  reception  failures,  which  leads  to  higher 
coordination  loss,  as  illustrated  in  Fig.  19(b).  This  results 
in  reduced  aggregate  throughput  seen  in  Fig.  19(c).  An 
interesting  observation  from  Fig.  19(b)  is  that  the  coordination 
loss  values  are  close  for  the  25  and  50  node  topologies. 
This  can  be  explained  as  follows.  Coordination  loss  occurs 
either  because  the  intended  receiver  is  transmitting  or  it 
is  beamformed  towards  a  different  direction  for  reception. 
Under  a  saturated  traffic  model,  each  node  always  has 
backlogged  packet  queues,  so  the  probability  that  a  node 
chooses  to  transmit  in  a  given  slot  is  independent  of  spatial 
density  of  nodes.  In  case  the  intended  receiver  chooses  not  to 
transmit,  the  probability  that  it  beamforms  towards  a  different 
direction  depends  on  the  probability  that  at  least  one  of  the 
other  neighbors  (in  addition  to  the  tagged  transmit  node) 
of  the  intended  receiver  chooses  to  transmit  to  that  node. 
The  increasing  node  density  has  two  opposing  effects  on 
this  probability:  the  probability  that  at  least  one  of  the  other 
neighbors  chooses  to  transmit  increases  with  node  density, 
but  the  probability  that  the  transmission  is  intended  for  the 
tagged  receiver  node  decreases  with  increasing  node  density 
(i.e.,  more  neighbors  for  each  node).  This  leads  to  similar 
coordination  losses  for  the  25  and  50  node  random  topologies. 

Note  that  the  following  two  points  are  important  while 
interpretting  the  directional  slotted  Aloha  results  in  light  of 
the  interference  analysis  presented  earlier:  (1)  our  interference 
analysis  finds  the  probability  of  packet  loss  due  to  interference, 
conditioned  on  the  assumption  that  the  tagged  transmitter  and 
the  receiver  are  beamformed  towards  each  other.  We  do  not 
consider  transmit-receive  coordination  loss  in  our  analysis. 
(2)  the  interference  loss  for  our  simulation  setting  is  expected 
to  be  lower  relative  to  that  obtained  for  the  flat-top  antenna 
and  linear  arrays  in  Section  IV-D  because  of  the  narrower 
azimuthal  half-power  beamwidth  of  the  patch  antenna 
array  (4°  versus  14.4°  and  8.5°:  given  the  relatively  smaller 
vertical  beamwidths  of  the  flat-top  antenna  and  also  the  flat-top 
elements  to  construct  the  linear  array,  the  required  azimuthal 
half-power  beamwidths  to  obtain  the  same  directivity 
are  larger).  In  other  words,  our  analytical  results  on  the 
interference  loss  are  more  conservative  than  the  simulations. 

The  low  interference  loss  results  for  the  slotted  Aloha 
example  indicate  that  the  pseudowired  abstraction  for  mm- 
wave  mesh  network  links  can  indeed  serve  as  an  appropriate 
first  order  approximation  for  MAC  design  in  highly  directional 
mm-wave  mesh  networks.  Therefore,  the  challenge  for  the 
network  designer  is  to  schedule  transmissions  in  the  face  of 
deafness  resulting  from  the  highly  directional  links,  while  ex¬ 
ploiting  the  reduced  interference  among  simultaneously  active 
mm-wave  links. 
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Fig.  19.  Simulation  results  for  directional  slotted  Aloha  protocol. 


VI.  Conclusions  and  Future  Work 


We  have  investigated  the  impact  of  the  unique  physical 
characteristics  of  mm-wave  links  on  the  statistics  of  spatial 
interference,  and  consequently  MAC  design  for  outdoor  mm- 
wave  mesh  networks.  Our  interference  analysis  framework 
enables  a  quantitative  evaluation  of  collision  loss  probability 
for  an  mm-wave  mesh  network  with  uncoordinated  transmis¬ 
sions,  as  a  function  of  the  antenna  patterns  and  spatial  density 
of  simultaneously  transmitting  nodes.  For  the  directivities 
typical  of  mm-wave  nodes  with  compact  form  factors,  the 
low  collision  loss  probabilities  indicate  that  the  pseudowired 
model  is  indeed  appropriate,  which  motivates  a  radically 
different  approach  to  MAC  design  relative  to  the  CSMA-based 
WiFi  mesh  networks  operating  in  the  low  frequency  bands. 
Rather  than  focusing  on  proactive  interference  management  as 
in  conventional  MAC  design,  we  must  now  devise  scheduling 
mechanisms  that  address  deafness. 

We  note  that  preliminary  results  [8]  for  Memory-guided 
Directional  MAC  (MDMAC),  a  MAC  protocol  designed 
according  to  these  guidelines,  are  promising.  However,  there 
are  a  number  of  important  topics  for  future  research  in 
the  design  of  mm-wave  mesh  networks.  It  is  important  to 
devise  automated  network  discovery  and  topology  update 
mechanisms,  closely  coupled  with  efficient  mechanisms  for 
transmit  and  receive  beamforming.  Another  key  issue  is  that  of 
packet-level  time  synchronization.  For  example,  the  MDMAC 
protocol  in  [8]  uses  learning  and  predictability  to  combat 
deafness,  but  is  predicated  on  network-wide  time  slotting. 
Lightweight  mechanisms  need  to  be  developed  for  acquiring 
and  maintaining  such  time  slotting  in  highly  directional 
networks.  We  also  need  more  detailed  propagation  models: 
while  rooftop-to-rooftop  links  may  be  well  approximated  as 
idealized  line-of-sight  links,  links  in  a  lamppost  based  network 
may  exhibit  severe  fading  even  with  the  sparse  multipath 
corresponding  to  a  highly  directional  link  [37].  This  has  an 
impact  both  on  physical  layer  design  (e.g.,  providing  antenna 
diversity)  and  network  architecture  (e.g.,  providing  route  diver¬ 
sity).  Finally,  design  for  indoor  mm-wave  mesh  networks  must 
address  additional  challenges  of  blockage  and  co-existence. 


Fig.  20.  Flat-top  antenna  pattern  in  polar  coordinates  (6,  cf))  around  the 
antenna  beam  axis. 


Appendix 


Antenna  directivity  computations:  We  use  (2)  in  order 
to  calculate  the  exact  antenna  directivity  from  the  antenna 
pattern.  We  first  need  to  compute  the  beam  solid  angle 
f l  =  fg_0  fJ=_7T  Gn{9,  </>)  sin  9  dd  d<j>  of  the  antenna,  where 
Gn(9,cj))  is  the  normalized  power  pattern. 

Flat-top  antenna:  Looking  at  the  3-dimensional  pattern  for 
a  fiat-top  antenna  with  sector-width  w,  in  polar  coordinates 
(9,  </>)  around  the  antenna  beam  axis  (see  Fig.  20),  we  have 


Gn{9A) 


1,  —7 r  <  <f>  <  tt,  0  <  9  <  |r, 
0,  otherwise. 


(34) 


Therefore, 


S2  = 


/  0—0  J  (f)——n 

C  j' 

1 0—0  J  tj)=—n 


Gn  ( 9 ,  </>)  sin  6  d9  d<j) 


sin  9  d9  d(j> 


=  2tt(1  -cos(-)). 


(35) 


Therefore,  directivity  D  of  the  flat-top  antenna  is  given  by 

4tt  2 


D  = 


27r(l  -  cos(f ))  l-cos(f)' 


(36) 
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(b)  Calculation  of  a 

Fig.  21.  Linear  array  of  flat-top  elements:  directivity  calculation. 


For  D  =  24  dBi,  we  obtain  w  =  14.4°,  which  is  the  required 
flat-top  antenna  sector-width  for  our  100m  mm-wave  mesh 
link  example. 

Linear  array  of  flat-top  antenna  elements:  For  a  simpler 
approach  to  calculation  of  the  beam  solid  angle  in  case  of  a 
uniform  linear  array  of  flat-top  elements  spaced  A/2  apart,  we 
re-orient  the  antenna  pattern  with  the  beam  axis  along  the  x- 
axis  (i.e.,  the  antenna  elements  in  the  array  are  placed  along  the 
y-axis),  as  shown  in  Fig.  21(a).  The  normalized  3-dimensional 
antenna  pattern  for  a  linear  array  of  N  flat-top  elements,  each 
with  a  beam-width  w,  in  this  coordinate  system  is  given  by 


Gn(M) 


{1  sin2  (^7r  cos  0) 
N*  sin2  (^7r  cos  0)  ’ 

0, 


7T  W  <- 

2  2  — 

otherwise. 


<f 


2  ’ 


(37) 


We  now  derive  ffO )  as  a  function  of  the  polar  angle  6.  Con¬ 
sider  the  circle  formed  by  the  intersection  of  the  3-dimensional 
flat-top  antenna  pattern  with  the  surface  of  the  (normalized) 
unit  sphere  centered  at  the  antenna  (Fig.  21).  From  Fig.  21(a), 
for  a  given  8 ,  the  limits  of  integration  for  <f>  are  given  by  ±cf>i 
where  <j>i  =  tan_1(|),  and  \ft\  <  |.  We  have  b  =  cos 
and  from  Fig.  21(b),  a(9)  =  (sin2(^)  —  cos 20)s.  Therefore, 

4>i{9)  =  tan-1  ^ - ^cob(^) — ~)-  We  now  compute  12  for 


a  linear  array  of  flat-top  elements  as  follows. 


J*>=  f- 


Consider  an  example  antenna  array  design  problem:  how  many 
flat-top  antenna  elements,  each  of  sector-width  20°,  spaced 
half  wavelength  (A/2)  apart  to  form  a  uniform  linear  array,  are 
needed  to  attain  directivity  D  =  24  dBi?  Using  the  calculations 
outlined  above,  we  find  that  a  linear  array  of  12  flat-top 
elements  provides  the  required  directivity.  We  use  this  linear 
array  for  the  100m  mm-wave  link  examples  in  the  paper.  Note 
that  we  can  also  derive  the  directivity  of  a  flat-top  antenna 
element  in  a  similar  fashion,  with  the  antenna  beam  axis  along 
the  x-axis  (i.e.,  Gn(0,0)  =  1  for  the  same  (9,cf>)  limits  as  in 
(37)). 
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